Accurate evolutions of inspiralling and magnetized neutron-stars: equal-mass binaries 
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By performing new, long and numerically accurate general-relativistic simulations of magnetized, equal-mass 
neutron-star binaries, we investigate the role that realistic magnetic fields may have in the evolution of these 
systems. In particular, we study the evolution of the magnetic fields and show that they can influence the survival 
of the hypermassive neutron star produced at the merger by accelerating its collapse to a black hole. We also 
provide evidence that, even if purely poloidal initially, the magnetic fields produced in the tori surrounding the 
black hole have toroidal and poloidal components of equivalent strength. When estimating the possibility that 
magnetic fields could have an impact on the gravitational-wave signals emitted by these systems either during 
the inspiral or after the merger, we conclude that for realistic magnetic-field strengths B < 10 12 G such effects 
could be detected, but only marginally, by detectors such as advanced LIGO or advanced Virgo. However, 
magnetically induced modifications could become detectable in the case of small-mass binaries and with the 
development of gravitational-wave detectors, such as the Einstein Telescope, with much higher sensitivities at 
frequencies larger than w 2 kHz. 

PACS numbers: 04.30.Db, 04.40.Dg, 04.70.Bw, 95.30.Qd, 97.60Jd 



I. INTRODUCTION 

The use of improved and more accurate numerical tech- 
niques, together with access to larger computational infras- 
tructures, has brought the simulation of binary neutron-star 
(BNS) systems to an unprecedented level of maturity. A num- 
ber of groups have reported on calculations of BNSs with dif- 
ferent levels of approximation, for equal- and unequal-mass 
systems, with and without magnetic fields (see, e.g., (T-8 1 for 
some of the most recent works). Besides the obvious impli- 
cations that these systems have in our understanding of the 
origin of short 7-ray bursts (GRBs), whose short rise times 
suggest that their central sources have to be highly relativis- 
tic objects [9|, BNS systems are expected to produce sig- 
nals of amplitude large enough to be relevant for Earth-based 
gravitational-wave (GW) detectors and to be sufficiently fre- 
quent sources to be detectable over the timescale in which the 
detectors are operative. The current estimate for the detection 
rate relative to the first-generation interferometric detectors is 
approximately 1 event per 40 — 300 years, increasing to an 
encouraging 10 — 100 events per year for the advanced detec- 
tors iflOl . 

The detection of gravitational waves from neutron-star (NS) 
binaries will also provide a wide variety of physical informa- 
tion on the component stars ifTTl . This includes their mass, 
spin, and radius, which would in turn provide vital clues on the 
governing equation of state (EOS), and, possibly, their mag- 
netic field. However, for this information to be extracted it 
is essential that accurate and long-term simulations are car- 
ried out, which span the interval ranging from the early in- 
spiral to the decaying tail of the late ringing of the formed 
black hole (BH). This is indeed the goal of this work, where 
we focus on whether or not present and future GW detectors 
will be able to determine the level of magnetization of NSs. 



This is not an academic question, as we know that NSs have 
very large magnetic fields, and it is indeed via the magnetic- 
dipolar losses that the vast majority of NSs are routinely de- 
tected as pulsars lfT2l . Yet, determining what the effects of 
magnetic fields are on the inspiral and merger of BNSs is a re- 
markably difficult task, requiring the solution of the Einstein 
equations together with those of general-relativistic magneto- 
hydrodynamics (GRMHD). So far, only three GRMHD simu- 
lations of inspiralling BNSs have been reported J3]|5]|6l, and 
while Refs. J3] |5) considered magnetic fields that are astro- 
physically unrealistic 1 |[T3llT4l . only the work in Ref. [6| has 
studied magnetic fields of the order of « 10 12 G, which are 
probably the strongest to be expected for NSs near the merger. 
Ultralarge magnetic fields are, however, not entirely uninter- 
esting from a general-relativistic point of view. Indeed, as 
discussed in [6], the magnetic tension associated with these 
extremely large magnetic fields can be so strong to reduce the 
stellar tidal deformations during the inspiral and hence to lead 
to a slightly delayed time of merger. 

Here we present a more extended analysis than the one 
given in [6 1 and report on a systematic investigation of equal- 
mass BNSs systems through long-term simulations using the 
highest resolutions to date. The calculations cover a range of 
magnetic fields from B ss 10 8 G up to B w 10 12 G, and two 
different masses to distinguish the phenomenology of those 
binaries that lead to a prompt collapse from those that lead 
instead to a delayed one (see the discussion in 0). Overall, 
we find that magnetic fields are amplified during the merger, 



We note that although NSs with magnetic fields as large as 10 are widely 
expected to be behind the phenomenology associated with magnetars, it is 
unrealistic to expect that the old NSs comprising the binary have magnetic 
fields that are so large. 
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when the turbulent motions, triggered during the merger by 
the Kelvin-Helmholtz instability, curl magnetic field lines pro- 
ducing a strong toroidal component that reaches a strength 
comparable to the poloidal one. The toroidal field main- 
tains a value comparable or larger than the poloidal one dur- 
ing the subsequent evolution of the hypermassive neutron star 
(HMNS) formed after the merger. The stability of the latter, 
however, is influenced by the strength of the poloidal field, 
which can transport the angular momentum outwards and trig- 
ger the collapse of the HMNS to a BH. Furthermore, equipar- 
tition among the poloidal and toroidal magnetic field compo- 
nents has been measured during the first 5 ms after the col- 
lapse of the HMNS, when the system consists of a rotating 
BH surrounded by a massive, high-density torus. 

We have also analyzed in detail the GW signal emitted 
by these systems and found that for the timescales consid- 
ered here, the overlaps in the GWs between a nonmagne- 
tized binary and a magnetized one are always above what 
detectors such as Advanced LIGO (advLIGO) or Advanced 
Virgo (advVirgo) can distinguish. Hence, it is very unlikely 
that present detectors will be able to measure the presence of 
magnetic fields. However, for sufficiently small-mass bina- 
ries, whose corresponding HMNS could survive for up to a 
fraction of a second (see the Appendix of [8 |), the dephasing 
induced by the presence of magnetic fields could be measur- 
able, especially by those detectors, such as the Einstein Tele- 
scope [ 15], that have higher sensitivities at frequencies larger 
than ss 2 kHz. 

The paper is organized as follows. In Sec. [II] we first 
summarize the formalism we adopt for the numerical solu- 
tion of the Einstein and of the GRMHD equations; we then 
describe briefly the numerical methods we implemented in 
the Whisky code lfl6HT8l . we outline our mesh-refined grid 
setup, and we finally describe the quasi-equilibrium initial 



data we use. In Sec. Ill we describe the dynamics of the dif- 



ferent models by studying both the evolution of the matter and 
of the magnetic field. In Sec. [iVjwe instead describe the GWs 
emitted by these systems and we estimate the possibility to 
detect magnetic field effects on those signals, while in Sec.[V] 
we summarize our main results. 

Here we use a spacelike signature (—,+,+,+) and a sys- 
tem of units in which c = G = M© = 1 (unless explicitly 
shown otherwise for convenience). 



II. MATHEMATICAL AND NUMERICAL SETUP 

Most of the details on the mathematical and numerical 
setup used for producing the results presented here are 
discussed in depth in |6, 18— 20|. In what follows, we limit 
ourselves to a brief overview and we describe in more 
details only the main differences with respect to our previous 
simulations. 

A. Einstein and Magnetohydrodynamics equations 



providing the solution of a conformal traceless formulation 
of the Einstein equations lfl9l . The GRMHD equations 
were instead solved using the Whisky code IT6HT81 , which 
adopts a flux-conservative formulation of the equations as pre- 
sented in ETTl and high-resolution shock-capturing schemes 
(HRSC). The Whisky code implements several reconstruc- 
tion methods, such as Total-Variation-Diminishing (TVD) 
methods, Essentially-Non-Oscillatory (ENO) methods [22] 
and the Piecewise Parabolic Method (PPM) J23). As al- 
ready discussed in [6| the use of reconstruction schemes of 
order high enough is fundamental for the accurate evolution 
of these systems and in particular for assessing the impact of 
the magnetic fields. Therefore all the results presented here 
have been computed using the PPM reconstruction, while the 
Harten-Lax-van Leer-Einfeldt (HLLE) approximate Riemann 
solver [24] has been used to compute the fluxes. 

In order to guarantee the divergence-free character of the 
MHD equations we have employed the flux-CD approach de- 
scribed in |25|, but with one substantial difference, namely, 
that we use as an evolution variable the vector potential in- 
stead of the magnetic field. In other words, by using an ex- 
pression similar to equation (31) of [25 ], we compute the elec- 
tric field at the center of each numerical cell by interpolating 
the fluxes computed at the interfaces of the cell and then use it 
to evolve directly the vector potential. We recall that in ideal 
MHD a relation exists between the fluxes of the magnetic field 



B and the value of the electric field E 



B l = ^B l 



-v X B, where 



(1) 

(2) 



and where 7 is the determinant of the 3-metric, v l is the 3- 
velocity of the fluid as measured by an Eulerian observer, a 
the lapse, and f3 l the shift vector. In particular, the following 
relations hold in Cartesian coordinates 



E x = F Z {B V ) = -F y (B z 



E,, 



-F Z (B X ) 
E z = F V (B X ) = 



F X {B Z ) , 
-F X (B V ) , 



with 



F l (B ] ) =i l B ] -VB 1 



(3) 
(4) 
(5) 



(6) 



The evolution equations for the vector potential A and for the 
magnetic field B can then be written as 



d t A = -E , 
B = V x A 



(7) 
(8) 



The evolution of the spacetime was obtained using the 
Ccatie code, a three-dimensional finite-differencing code 



Equation (|7]i is solved at the center of each cell (i, j, k), where 
the electric field is given by 
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E x { Xi , yj ,z k ) = -(-F y (B 2 ) (jJ+1/2 , fe) -^(S z ) (lJ „ lAfc) +^(^) ( ^ , (9) 

Ey(%i,yj,Zk) = ^ (F x {B z )(i+i/2^,k) + F x {B z )(i-\/2,j,k) - F z {B x )(i t j,k+\/2) - F z (B x )(i,j,k-i/2)j > (10) 
E z ( Xi , yj ,z k ) = J(-F*(^) (i+1 / 2J ,fc)-F*(^ , (ID 



F % (Bi ) being the numerical flux computed at the interface of 
the cell. 

Since the magnetic field is computed from the curl of the 
vector potential using the same differential operator used to 
compute its divergence {i.e., a central-difference scheme), 
its divergence free character is guaranteed at essentially ma- 
chine precision at all times, also when using adaptive mesh- 
refinement (AMR). We note that a similar approach has been 
recently implemented also in another code 11261 and, in anal- 
ogy with |26l . we add a Kreiss-Oliger type of dissipation E71 
to the evolution equation of the vector potential in order to 
avoid the possible formation of spurious post-shock oscilla- 
tions in the magnetic-field evolution. It has indeed been shown 
by 1 28 1 that applying TVD operators to the vector potential 
does not guarantee automatically the TVD character of the 
magnetic field, leading to possible post-shock oscillations in 
the latter. The code has been validated against a series of 
tests in special relativity [29] and in full general relativity 
(seed]). 

The system of GRMHD equations is closed by an EOS and, 
as discussed in detail in [2], the choice of the EOS plays a fun- 
damental role in the post-merger dynamics and significantly 
influences the survival time against gravitational collapse of 
the HMNS produced by the merger. 

As already done in (6), also in this paper we have employed 
the commonly used "ideal-fluid" EOS, in which the pressure 
p is expressed as p = pe{T — 1), where p is the rest-mass 
density, e is the specific internal energy and T is the adiabatic 
exponent. Such an EOS, while simple, provides a reasonable 
approximation and we expect that the use of realistic EOSs 
would not change the main results of this work. 

B. Adaptive Mesh Refinements 

Both the Einstein and the GRMHD equations are solved 
using the vertex-centered AMR approach provided by the 
Carpet driver [30]. Our rather basic form of AMR consists 
in centering the highest-resolution level around the peak in 
the rest-mass density of each star and in moving the "boxes" 
following the position of this maximum as the stars orbit. The 
boxes are evolved as a single refinement level when they over- 
lap. 

The results presented below refer to simulations performed 
using 6 levels of mesh refinement with the finest level having 
a resolution of h = 0.1500 M ~ 221m. The grid structure 
is such that the size of the finest grids is 24 M & ~ 35.4 km, 
while a single refinement level covers the region between a 



distance r = 164 M ~ 242.2 km and r = 254.4 Mq ~ 
375.7 km from the center of the domain. This region is the 
one in which our gravitational-wave extraction is carried out, 
with a resolution of h = 4.8 Mq ~ 7.1 km (as a comparison, 
the gravitational wavelength is about 100 km and thus well- 
resolved on this grid). In addition, a set of refined but fixed 
grids is set up at the center of the computational domain so as 
to better capture the details of the Kelvin-Helmholtz instabil- 
ity (cf. (5]). Moreover, after the merger, at about 8.5ms, we 
enlarge the central grid that is formed by the merging of the 
two initial boxes. We do this in order to cover a cubical region 
with a side of about 88.6 km and so better resolve not only 
the whole HMNS, but also the BH-torus system which is pro- 
duced by the collapse of the HMNS. For all the simulations 
reported here we have used a reflection-symmetry condition 
across the z = plane and a 7r-symmetry condition across the 
x = plane 2 . At the outer boundary we instead used simple 
zeroth-order extrapolation on the MHD variables (in practice, 
we just copy the value of the MHD quantities from the outer- 
most evolved point in each direction to the points of the outer 
boundary in that direction). Also note that a very little amount 
of matter and magnetic fields reaches the outer boundary, so 
the effect of the outer-boundary conditions on the MHD and 
hydrodynamical variables is negligible. 

The timestep on each grid is set by the Courant condition 
(expressed in terms of the speed of light) and so by the spa- 
tial grid resolution for that level; the Courant coefficient is 
set to be 0.35 on all refinement levels. The time evolution is 
carried out using 4th-order-accurate Runge-Kutta integration 
algorithm. Boundary data for finer grids are calculated with 
spatial prolongation operators employing 3rd-order polyno- 
mials for the matter variables and 5th-order polynomials for 
the spacetime variables. The prolongation in time employs 
2nd-order polynomials and this ensures a significant memory 
saving, requiring only three timelevels to be stored, with little 
loss of accuracy due to the long dynamical timescale relative 
to the typical grid timestep. 

The grid setup used here is therefore quite different from 
the one adopted in our previous work on magnetized NS bi- 
naries [6], where we used fixed mesh-refinement in order to 
reduce the violation (generated by the interpolation in the 
buffer zones) of the divergence-free constraint of the mag- 
netic field. Our current implementation, based of the evolu- 



2 Stated differently, we evolve only the region {x > 0, z > 0} applying a 
180° -rotational-symmetry boundary condition across the plane at x = 0. 
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TABLE I: Properties of the eight equal-mass binaries considered: proper separation between the stellar centers d/M ADM ; baryon mass Mb 
of each star; total ADM mass A/ ADM ; angular momentum J; initial orbital angular velocity Qo~, mean coordinate radius r e along the line 
connecting the two stars; ratio of the polar to the equatorial coordinate radii r p /r e ; maximum rest-mass density p m a X ; maximum initial 
magnetic field Bo, where * is 8, 10 or 12. Note that M ADM and J are reported as measured on the finite-difference grid. 



Binary 


d/M ADU 


M b (M ) 


M AOM (M Q ) 


J (gcm 2 /s) 


flo (rad/ms) 


r e (km) 


r p /r e 


Pmax (gm/cm 3 ) 


Bo (G) 


Ml .45-B* 


14.4 


1.445 


2.680 


6.5084 x 10 49 


1.78 


15.0 ±0.3 


0.899 


4.58 x 10 14 


or 1.97 x 10* 


Ml . 62-B* 


13.3 


1.625 


2.981 


7.7806 x 10 49 


1.85 


13.6 ±0.3 


0.931 


5.91 x 10 14 


or 1.97 x 10* 



tion of the vector potential, does not produce any violation of 
the divergence-free condition of the magnetic field, since it in- 
terpolates the vector potential instead of the magnetic field in 
the buffer zones. Moreover, since the vector potential is stored 
at the center of the cell, it is possible to use without modifica- 
tion the prolongation and restriction operators currently avail- 
able in the Carpet driver. This makes it possible to use the 
moving-grid setup that has been utilized with success in our 
previous general-relativistic hydrodynamics simulations. 

C. Initial data 

The initial data are the same as those used in (2] |6j. 
They were produced by Taniguchi and Gourgoulhon OTI 
with the multi -domain spectral-method code LORENE (32). 
The initial solutions for the binaries are obtained assuming a 
quasi-circular orbit, an irrotational fluid-velocity field, and a 
conformally-flat spatial metric. The matter is modeled using 
a polytropic EOS p = Kp T with K = 123.6 and T = 2, 
in which case the maximum gravitational mass is M ADM ~ 
1.82 Mq for a nonrotating star and A/ ADM ~ 2.09 M for 
a uniformly rotating one. Since no self-consistent solution 
is available for magnetized binaries yet, a poloidal magnetic 
field is added a-posteriori using the vector potential 

A 4> = zu 2 A b max(p~p cut ,0) ns , (12) 

where w = y^x 2 + y 2 , A b > parameterizes the strength 
of the magnetic field, p cut defines where in the NS the mag- 
netic field goes to zero, and n B determines the degree of differ- 
entiability of the potential. The components of the magnetic 
field are then computed by taking the curl of the Cartesian 
components of Eq. ( p"2] i to enforce that the divergence of the 
magnetic field is zero at machine precision. Here we have set 
Pcut = 0.04max(P), and n s = 2 to enforce that both the 
magnetic field and its first derivative are zero at p = p cu t- In 
Ref. OJ the magnetic field was built with an expression equiv- 
alent to ( p"2| ), but with p cut set to the pressure in the atmo- 
sphere, and in Ref. [5 1 the expression used is slightly different 
and P cut is set to be 4% - 0.1% of max(p); in both Refs. [3 1 
and n s = 1. 

Table|I]lists some of the properties of the eight equal-mass 
binaries considered here. More specifically, we have consid- 
ered two classes of binaries differing in the initial masses, 
i.e., binaries Ml . 4 5-B*, and binaries Ml . 62-B*. For each 
of these classes we have considered four different magne- 
tizations (indicated by the asterisk) so that, for instance, 



Ml . 4 5-B12 is a low-mass binary with a maximum initial 
magnetic field B = 1.97 x 10 12 G. Note that the binaries 
with zero magnetic fields are the same as those evolved in 
Ref. 0. 



D. Gravitational- Wave Extraction 

Details about the algorithms implemented in the code to ex- 
tract the GW signal can be found in (2)- Here we just remind 
the reader that we compute the waveforms using two different 
methods. The first one is based on the Newman-Penrose for- 
malism and computes the Weyl scalar \&4. The gravitational- 
wave polarization amplitudes h + and h x are then related to 
*5>4 by simple time integrals 03l 

h + -ih x =V 4 , (13) 

where the double overdot stands for the second-order time 
derivative. 

The second method is instead based on the measurements of 
the nonspherical gauge-invariant perturbations of a Schwarz- 
schild BH (see refs. 0414361 for some applications of this 
method to Cartesian-coordinate grids). In practice, a set of 
"observers" is placed on 2-spheres of fixed radius where we 
extract the gauge-invariant, odd-parity (or axial) current mul- 
tipoles Q^ m and even-parity (or polar) mass multipoles Q\ m 
of the metric perturbation 0711381 . The Qt m an d Qi m vari- 
ables are related to h+ and h x as [39 1 

h+ ih x = -L J2 (QL i QL Wj -2Y lm . 

(14) 

Here _2Y tm are the s = — 2 spin- weighted spherical harmon- 
ics and (£, m) are the indices of the angular decomposition. 

Since the two methods have been shown to give wave- 
forms that are identical up to the truncation error, we will 
here use h + computed only with the gauge-invariant quan- 
tities and we will focus only on the I = 2, to = 2 mode 
since the others have amplitudes which are negligible com- 
pared to this. All the waveforms have been extracted at a ra- 
dius 7i SO = 2OOM 300 km. We also ignored the contri- 
bution from the spherical harmonics since they depend on the 
direction of the source with respect to the detector and con- 
tribute as a multiplication factor of order 1; thus they do not 
modify the results presented here. 




FIG. 1: Snapshots at representative times of the evolution of the high-mass binary with initial maximum magnetic field of 10 10 G, 
i.e., Ml . 62-B10. Shown with two different color-code maps are the rest-mass density p (red-yellow) and the magnetic field \B\ (blue- 
green-yellow- white). To better visualize the inner structure we plot only the values on z < 0. In order to show the two scalar quantities at the 
same time, they are shown on either side of a fictitious screen (p on the left and |Bj on the right). The first four panels refer respectively to 
the binary at the initial separation of 45 km, to the binary after two orbits, to the merger and to the bar-deformed HMNS. The last two panels, 
instead, refer respectively to when the BH has just been formed and to a subsequent stage of the quasi-stationary evolution of the BH-torus 
system. The grey spheroidal surface in the center represents the location of the apparent horizon. 
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E. Accuracy of the Results 

A reliable assessment of the truncation error is essential to 
draw robust conclusions on the results of numerical simula- 
tions. Following a procedure discussed in detail in Ref. PI , 
also here we have carried out a systematic measurement of 
the accuracy and convergence properties of our simulations, 
and deduced a corresponding "error-budget". The main con- 
clusions are very similar to those drawn in Ref. [4|, which 
for compactness only we briefly recall here. More specif- 
ically, we showed that with typical (finest) resolutions of 
h ~ 0.12 M© — 0.19 Mq, the results show the expected con- 
vergence rate of 1.8 during the inspiral phase, which how- 
ever drops to 1.2 at the merger and during the evolution of 
the HMNS. This deterioration of the convergence rate is due 
mostly to the strong shocks which form during the merger and 
which HRSC schemes can reproduce at lst-order only. Fur- 
thermore, physical quantities, such as the rest-mass, are con- 
served with a relative error of < 10~ 6 , while the energy and 
the angular momentum are conserved to < 1% after taking 
into account the parts lost to radiation. Finally, the expected 
agreement in both phase and amplitude is found in the wave- 
forms extracted from different detectors within the same sim- 
ulation or from the same detector but at different resolutions. 
Such waveforms have been found to be also convergent at a 
rate of 1.8 (see 0J for details). Finally, for all the simulations 
reported here the violation of the Hamiltonian constraint has 
an L2-norm which is < 10- 4 /A/| DM for the high -mass bina- 
ries and < 10 _5 /Af^ DM for the low-mass ones, for which no 
BH is formed. 



III. BINARY DYNAMICS 

As mentioned above, in order to highlight some of the 
most salient aspects of the binary dynamics it will be suffi- 
cient to consider two main classes of initial configurations: 
Ml . 62 -B* and 1 . 45-B*. These models differ only in the 
mass, the first being composed of stars each having a rest mass 
of 1.625 Mq (which we refer to as the "high-mass binaries"), 
the second of stars of rest mass 1.445 Af (which we refer to 
as the "low-mass binaries"). The use of these two classes 
is useful to distinguish the phenomenology of binaries whose 
merger leads to a prompt collapse of the HMNS from those 
where the HMNS can instead survive for several tens of mil- 
liseconds and up to a fraction of a second (see the discussion 
in J2))- We also note that in the case of the unmagnetized 
models, the dynamics is the same as the ones described in |2), 
to which we refer the interested reader for a more detailed 
description of the evolution of the matter and of the hydrody- 
namical instabilities such as the Kelvin-Helmholtz instability. 

A synthetic overview of the dynamics is summarized in 
Fig. [T] which shows snapshots at representative times of the 
evolution of the high-mass binary with an initial maximum 
magnetic field of 10 10 G, i.e., Ml . 62-B10. Shown with two 
different color-code maps are the rest-mass density p (red- 
yellow) and the magnetic field \B\ (blue-green-white). To 
better visualize the inner structure we plot only the values 



on z < 0. In order to show the two scalar quantities at the 
same time, they are shown on either side of a fictitious screen 
(p on the left and \B\ on the right). The first four panels re- 
fer respectively to the binary at the initial separation of 45 
km (t — ms), to the binary after two orbits (t = 4.7 ms), 
to the merger {t — 9.2 ms) and to the bar-deformed HMNS 
(t = 10.6 ms). The last two panels, instead, refer respectively 
to when the BH has just been formed (t = 12.6 ms) and to a 
subsequent stage of the quasi-stationary evolution of the BH- 
torus system (t = 15.2 ms). 

With this overall qualitative behavior of the binary in mind, 
we will next consider a more quantitative discussion of the 
evolution of the magnetic fields and we will only briefly sum- 
marize the dynamics of the matter. In doing this we will 
present in Figs.|2]and|3]the evolution of both the high and low- 
mass binaries to aid the comparison between the two classes 
of models. 



A. High-mass binaries 

We start by considering the evolution of the high-mass bi- 
naries Ml . 62 -B*, some of which were already considered 
in (5), where it was shown that initial magnetic fields lower 
than 10 14 G do not affect the dynamics in the inspiral phase. 
Overall, given the initial coordinate separation of 45 km, all 
binaries inspiral for approximately 3 orbits before merging at 
( « 8.2 ms. There are different ways to measure the time of 
the merger and the one we adopt here consists in looking at the 
first peak in the evolution of |^4|. This time corresponds ap- 
proximately to when the two stellar cores merge and we note 
that the external layers of the stars enter into contact about 
2 ms earlier. In the top panel of the left column of Fig. [2] we 
show the evolution of the maximum of the rest-mass density 
Pmax normalized to its initial value. It is particularly clear 
from the evolution of p max that all the models merge at the 
same time (e.g. see the first minimum in the evolution), while 
the post-merger dynamics are quite different. All the mod- 
els form an HMNS that survives a few milliseconds before 
collapsing to a Kerr BH, but its survival time varies consid- 
erably, as well as the number of oscillations in the evolution 
of the density before the rapid exponential increase in corre- 
spondence with the collapse. A discussion about this will be 



presented in Sec. Ill B 



The middle panel of the left column of Fig. [2] shows instead 
the maximum of the absolute value of the divergence of the 
magnetic field. To the best of our knowledge this is the first 
time that the evolution of the divergence of the magnetic field 
is shown in a GRMHD simulation of BNSs. Because this is 
a fundamental quantity to evaluate the quality of a numerical 
calculation, we encourage other authors to present it systemat- 
ically as well. As expected on mathematical grounds, the im- 



plementation of the GRMHD equations discussed in Sect. II A 



is such that the divergence of the magnetic field is essentially 
at machine precision at all times. It is important to stress that 
such a small violation would not be possible with the cell- 
centered AMR algorithm provided by the Carpet code un- 
less the vector potential is used as an evolved variable. 
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FIG. 2: Evolution of the maximum of the rest-mass density p normalized to its initial value (top row), of the maximum of the absolute value 
of the divergence of B (middle row), and of the maximum of the magnetic field strength \B\ (bottom row). The left and right columns refer to 
the high-mass and low-mass binaries, respectively. Note that in the case of the high-mass models (left column), the values of |B ma x| after BH 
formation refer, for the large majority of the time, to matter outside the apparent horizon and in the torus. 
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Finally, the bottom panel of the left column of Fig. [2]shows 
that the magnetic field grows mostly at the time of the merger 
and reaches values which are about one order of magnitude 
higher, before the collapse to BH. We note that in the case 
of the high-mass models (left column) the values of |-B max | 
after BH formation refer, for the large majority of the time, 
to matter outside the apparent horizon and in the torus. This 
is because the steep gradients of the matter variables inside 
the apparent horizon are under-resolved as a result of the grid 
stretching and dissipated on a timescale which is of the order 
of a fraction of a ms. Hence, with a few possible exceptions, 
the data in the plots refers statistically to the matter outside 
the apparent horizon. 

The growth of the magnetic field at the merger is made more 
clear in the different panels contained in the left column of 
Fig. [3] where we concentrate in particular on the evolution of 
the maxima of the total magnetic field (black solid line) and 
of its toroidal (red dot-dashed line) and poloidal (blue long- 
dashed line) components. As already shown in [2|, Kelvin- 
Helmholtz instability develops during the merger, when the 
external layers of the two NSs enter into contact, i.e., roughly 
2 ms before the time of the merger, which is indicated in those 
panels with the first vertical dotted line. This purely hydrody- 
namical instability leads to the formation of vortices that can 
curl magnetic field lines that were initially purely poloidal and 
produce toroidal components. As it is evident from the panels 
in Fig. [3] a strong toroidal component is indeed formed in all 
cases and it reaches values that are comparable or larger than 
the poloidal component, but its energy is not in equipartition 
with the kinetic energy in the layer. Despite the exponential 



growth caused by the Kelvin-Helmholtz instability, the overall 
amplification of the magnetic field is of an order of magnitude 
at most, with a growth rate dB/dt ~2x 10 12 (G/ms) in the 
case of model Ml . 62-B12. This is in contrast with what was 
reported by 14011 . where an amplification of several orders of 
magnitude in the magnetic field of the HMNS was observed, 
with a growth rate dB/dt ~ 2 x 10 15 (G/ms) for a model 
similar to Ml . 62-B12. 

It is presently unclear what the origin of this discrepancy 
is. It is possible that this is due to the use of very dif- 
ferent numerical techniques, namely smooth-particle hydro- 
dynamics and HRSC methods. It is also possible that al- 
though we have used the largest resolutions employed so 
far in simulations of magnetized BNSs, such resolutions are 
not yet sufficient to properly resolve the nonlinear develop- 
ment of the instability. Studies of the effect of these insta- 
bilities and of the consequent amplification of the magnetic 
fields have recently been performed with local simulations 
on simpler backgrounds BT1 l42l . These studies have indeed 
shown that in order to achieve convergence in the vortex re- 
gion it is necessary to use resolutions that are much higher 
than those currently affordable in BNS simulations. On the 
other hand, by performing simulations with different resolu- 
tions for model Ml . 62 -B 12 we did not observe any sensi- 
ble difference in the amplification of the magnetic field and 
indeed the magnetic field evolution is certainly consistent if 
not convergent (see the discussion in |4| about why it is dif- 
ficult to determine the convergence order after the merger). 
This is shown in Fig.[4j where we report the evolution of the 
maximum magnetic field strength |-B| max for the high-mass 
model M1.62-B12 evolved with three different resolutions: 
h = 177 m (high resolution, black solid line), h = 221m 
(medium resolution, which is the standard resolution used in 
this article, blue short-dashed line) and h = 354 m (low res- 
olution, red long-dashed line). The curves have been shifted 
in time to account for the slightly different time of the merger. 
It is clear that doubling the resolution produces a difference 
in the amplification of less than a factor of about 2 (compare 
the red long-dashed line with the black solid line); the differ- 
ences become even smaller when comparing the medium and 
high resolution 3 . A similar consistency with resolution is not 
present in the simulations reported in |40l , where the differ- 
ences among the amplified magnetic fields seem to become 
even larger with increasing resolution 4 . 

Overall we believe that the main reason why the toroidal 
magnetic field in our simulations does not grow significantly 
at the merger is that the timescale over which the instabil- 
ity can develop is rather short. The shear layer between the 
two stars, in fact, survives only for about 1 ms, before being 



3 Figure [4] also shows a considerable increase in the magnetic field at the 
merger. However, this is not related to the Kelvin-Helmholtz instability, 
but rather to flux conservation which amplifies the magnetic field when the 
matter is compressed by the collision of the two stellar cores. 

4 Note that because we are here capturing a non-sustained turbulent flow, the 
variations of the magnetic field strength with resolution are not necessarily 
monotonic. 
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destroyed by the collision between the two stellar cores. In 
Refs. ll4Tll42l . the amplification of the magnetic field has been 
observed on timescales that are even shorter than this one, but 
only under very specific conditions and only for specific val- 
ues of the velocity at the shear layer. The differences between 
the condition under which the instability develops in our fully 
general-relativistic simulations and those used in these local 
simulations may explain the different results. Clearly, the best 
way to assess whether or not the development of the Kelvin- 
Helmholtz instability leads to a large or only to a moderate 
field amplification is to perform direct comparisons with other 
general-relativistic simulations of magnetized BNSs. Unfor- 
tunately, so far the only other reported evolution of the mag- 
netic field is the one in O, which is clearly not useful for an 
independent comparison. 

As a final remark, it is important to emphasize that the 
toroidal and poloidal components have comparable values 
also in the torus that is formed after the collapse to BH 
(cf. panels in the left column of Fig. [3]). Since most of the sim- 
ulations to date of magnetized accretion disks around BHs that 
model the central engine of short GRBs use initial conditions 
in which the magnetic field has only a poloidal component, it 
is of particular importance to remark that more realistic initial 
data should instead have a toroidal and a poloidal component 
of comparable magnitude. 



B. Delay of the Collapse 
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FIG. 5: Lifetime of the HMNS formed after the merger in the high- 
mass case as a function of the initial magnetic field. The error bar 
has been estimated from a set of simulations of unmagnetized binary 
NS mergers at three different resolutions; in particular, we have as- 
sumed that the magnetized runs have the same relative error on the 
delay time of the corresponding unmagnetized model. Indicated with 
a dashed line is the continuation of the delay times to ultra-high mag- 
netic fields of 10 17 G. 



It has been shown and discussed in a number of works that 
when the merger leads to a collapse, the time of survival of the 
HMNS depends on several factors, which include: the EOS, 
the efficiency in the redistribution of angular momentum, and 
the efficiency of the radiative transfer. Clearly, all of these 
influencing factors will act differently in highly-magnetized 
matter and hence the delay time t^, i.e., the time between 
the formation of the HMNS and its collapse to a BH, can be 
used to measure indirectly the magnetic fields of the progen- 
itor NSs. There are several different ways of defining tj, but 
a convenient and gauge-invariant one is to consider the delay 
time as the interval between the first and last peak in the evo- 
lution of 1 , which are always well-defined in the amplitude 
evolution, as these can be taken to correspond to the merger 
of the stellar cores and to the BH production. 

In Fig. [5] we show therefore the survival time of the HMNS 
as a function of the initial magnetic field strength, together 
with the error bar as estimated from a set of simulations of 
unmagnetized binary NS mergers at different resolutions (the 
delay time converges at first order, increasing with resolu- 
tion). It is clear from Fig. [5]that while models Ml . 62-BO and 
Ml . 62 -B8 have roughly the same post-merger dynamics and 
the same collapse time (see also the top left panel of Fig. [2]), it 
is also clear that models Ml . 62-B10 and Ml . 62-B12 col- 
lapse earlier than the unmagnetized one. To understand why 
this is the case, we recall that magnetic fields can affect the 
dynamics of the HMNS as first shown in axisymmetric evo- 
lutions of an isolated differentially rotating HMNS [43, 44 1. 



In essence, magnetic fields can, via magnetic tension 5 , redis- 
tribute the angular momentum, transporting it outwards and 
reducing the amount of differential rotation that is essential 
in supporting the HMNS against gravitational collapse (we 
recall that a HMNS has, by definition, a mass which cannot 
be sustained by the star if rotating uniformly). The ratio be- 
tween the magnetic tension and the pressure gradients scales 
like the ratio between the magnetic pressure and the gas pres- 
sure, and this ratio increases (although remaining less than 
one) after the merger because the magnetic fields are stronger 
and the HMNS is more extended and has smaller pressure gra- 
dients. As a result, magnetic fields can "accelerate" the col- 
lapse of these models, but only if they are sufficiently strong 
so that the magnetic tension can be comparable to or larger 
than the normal pressure gradients. Hence, the efficiency in 
angular-momentum redistribution will be proportional to the 
intensity of the (square of the) magnetic field and this explains 
why the delay time is essentially unchanged for small mag- 
netic fields, such as Bq < 10 s G. For larger values, how- 



5 We recall that in Newtonian ideal MHD the Lorentz force appearing in the 
equation for the conservation of momentum is given by 



1 

4np 



(V X B) X B 



1 

AfKp 



(B ■ V)B — V 



(15) 



where in the right-hand side the first term is the "magnetic tension" along 
the field lines and the second one is the (isotropic) "magnetic pressure". 
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FIG. 6: Evolution of the total rest-mass for the different high-mass 
binaries considered. Note that the sudden drop corresponds to when 
the apparent horizon is formed since we exclude the region inside it 
from the computation of the mass. Note also that the early collapse 
of Ml . 62-B10 leads to tori which are about a factor of two less 
massive. 



ever, the magnetic fields can influence the dynamics of the 
HMNS and decrease Td as shown by models Ml . 62 -BIO 
and Ml . 62-B12. 

Interestingly, the HMNS relative to the binary 
Ml . 62-B12 survives longer than the Ml . 62 -BIO one. This 
should not be entirely surprising since a very large magnetic 
field will also introduce a magnetic pressure [cf. eq. (fT5j)], 
which will provide an additional pressure support and thus 
either compensate or even dominate the angular-momentum 
redistribution. Indeed, when simulating a binary with an 
initial magnetic field of Bo ~ 10 17 G we have found that the 
delay time increases and is even larger than the one obtained 
in the absence of a magnetic field. This is not shown directly 
in Fig. [5] which has been restricted to realistic values of 
the magnetic field, but we have indicated with a dashed line 
the continuation of the delay times to ultra-high magnetic 
fields. Clearly, because of this tight correlation between the 
degree of magnetization of the NS matter and the delay of the 
time of the collapse, the measurement of the latter via a GW 
detection will allow to infer the former. 

The difference in the time of the collapse produces also 
small differences in the mass of the final BH and torus. This 
is shown in Fig. [6| which reports the evolution of the total 
rest-mass for the different high-mass binaries considered, and 
where the sudden drop corresponds to the formation of the ap- 
parent horizon (the matter inside the horizon is excluded from 
the computation of the baryon mass; see P31 for a discus- 
sion of the properties of the collapse with the gauge condi- 
tions used here). Similarly, in Table [il] we list the mass and 
spin of the BH formed at the end of the evolution, and the 
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FIG. 7: Comparison of the properties of the tori produced either by 
a magnetized binary (Ml . 62-B12, blue dashed line) or by a un- 
magnetized one (Ml . 62-BO, red solid line). Top panel: rest-mass 
density along the a;-axis at about 3 ms after the formation of the ap- 
parent horizon and which we truncate at 10 10 g/cm 3 . Bottom panel: 
Angular velocity at the same time as above; shown as reference with 
a dotted line is the Keplerian angular velocity S7kc P , which matches 
very well the outer parts of the torus. 



mass and radius of the torus. Since the models collapse at 
different times we have taken our measure at the end of the 
simulation (i.e., at t ~ 20 ms), when the accretion onto the 
BH is small and essentially stationary. In all cases the mass of 
the BH is M B h ~ 2.9M and the spin is a = J/M 2 « 0.8, 
but the mass of the torus drops from about 0.063 — 0.085 M Q 
to 0.033M© in the case of model Ml . 62-B10. This is prob- 
ably due to the fact that the magnetic field causes some matter 
to move outside the core region and that will become a BH; 
as a result, the longer the delay time, the larger the tori. We 
note that, even if small, these tori could still provide sufficient 
energy to power short GRBs. 

As a final remark we note that at least over the timescales 
considered here, the differences in the local dynamics of the 
torus matter between magnetized and unmagnetized binaries 
is very small. This is because the magnetic field is not yet 
strong enough to produce significant changes in the dynam- 
ics. A convincing example is shown in Fig. [7] which offers 
a comparison of the properties of the tori produced either by 
a magnetized binary (Ml . 62-B 12, blue dashed line) or by 
a unmagnetized one (Ml . 62-B0, red solid line). The top 
panel, in particular, shows the rest-mass density along the x- 
axis at about 3 ms after the formation of the apparent horizon 
and which we truncate at 10 10 g/cm 3 . Besides small differ- 
ences (the data refers to very different simulations), the den- 
sity profiles are very similar. An analogous conclusion can 
be drawn when looking at the bottom panel, which shows the 
angular velocity at the same time as above; also reported as 
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TABLE II: Columns 2 — 3 report the mass M and spin a of the 
BH, while column 4 shows the mass of the torus formed after the 
merger of the high-mass models, and column 5 the radius of the torus 
(computed as a mean of the position where the rest-mass density goes 
below 10 10 gcm~ 3 in the time interval between 19 and 20 ms). All 
the other quantities have been measured at t — 20 ms, when the 
accretion onto the BH is small and essentially stationary. 



Binary 


M [M©] 


a = J/M 2 


M tor [M ] 


r tor [km] 


Ml . 62-BO 


2.90 


0.80 


0.076 


105 ±13 


Ml . 62-B8 


2.89 


0.80 


0.085 


102 ± 16 


Ml . 62-B10 


2.94 


0.82 


0.033 


69 ±4 


Ml . 62-B12 


2.91 


0.81 


0.063 


94 ±4 



reference with a dotted line is the Keplerian angular velocity 
^Kep, which matches very well the outer parts of the torus. 



C. Low-mass binaries 

As already shown in |!6 | and anticipated in the previous 
Section, also in the low-mass case the presence of an initial 
magnetic field introduces no significant modification in the 
evolution of the binaries during the inspiral. To compare di- 
rectly with the behavior of the high-mass binaries, we show 
in the right column of Fig.|2]the evolution of the maximum of 
the rest-mass density normalized to its initial value, the max- 
imum of the absolute value of the divergence of the magnetic 
field and the maximum of the magnetic field. We recall that 
in Ref. [8] it was shown that the low-mass models take more 
than 100 ms to collapse to BH in the unmagnetized case, so 
the 20 ms of evolution of the present work are not sufficient to 
reach the collapse. 

In the top panel it is possible to appreciate that the evolu- 
tion of p max for the unmagnetized case and those of the mag- 
netized binaries are very similar, with only small differences 
in the frequency of the oscillations of the HMNS formed af- 
ter the merger. Overall, the presence of a magnetic field de- 
creases the oscillation frequency (cf. inset), probably because 
the additional magnetic tension counters the expansions of the 
bar-deformed HMNS. Moving over to the right bottom panel 
of Fig. [2] it is possible to note that also in the low-mass case all 
the magnetized models show an amplification of the magnetic 
field of about one order of magnitude and also in this case the 
divergence of the magnetic field is zero essentially at machine 
precision. As for the high-mass binaries, interesting point to 
note is that, on the timescale studied here, the magnetic field 
grows by about one order of magnitude soon after the merger 
(at t ps 8 ms), but then it saturates to a constant value. 

Additional information about the magnetic -field evolutions 
for the three different models are given in the panels in the 
right column of Fig. [3] Also in this case the toroidal compo- 
nent of the magnetic field (red dot-dashed line) is amplified 
exponentially because of the Kelvin-Helmholtz instability at 
the time of the merger of the external layers of the stars and 
it reaches the same value of the poloidal component. Both 
components have comparable values for the remaining dura- 



tion of the simulation and we expect that also in this case the 
collapse of the HMNS will produce a torus with a magnetic 
field configuration in which the toroidal and poloidal compo- 
nents have the same strength. This seems to be, at least for the 
equal-mass BNSs considered here, a universal characteristic 
of the tori that are formed from these systems. 



IV. GRAVITATIONAL-WAVE EMISSION 

A. High-mass binaries 

In Fig. [8] we show the GW signals emitted by the 4 high- 
mass binaries considered in this paper. The top left panel 
shows the unmagnetized case, the top right panel the model 
with an initial magnetic field of 10 8 G, the bottom left panel 
B w 10 10 G and the bottom right panel B ps 10 12 G. In the 
bottom right panel, together with M1.62-B12 (black solid 
line, which collapses at t ps 16 ms), we also show - to make 
the comparison clearer - the evolution of Ml. 62-BO (red 
dashed line, which collapses later). All the waveforms ex- 
hibit very similar features and, with the exception for the dif- 
ferent duration of the post-merger phase already discussed in 



Sec. Ill A they are almost indistinguishable from each other. 
Therefore, for all the models the signal is essentially com- 
posed of three parts: the inspiral (from t — r = ms to 
t — r ps 8 ms), the HMNS evolution (from t — r ps 8 ms 
to t — r ps 13 — 17 ms) and the ring-down of the final BH. 
The high-frequency oscillations in the post-merger phase are 
due to the cores of the two NSs that repeatedly bounce against 
each other until a sufficient amount of angular momentum is 
extracted via GWs emission or is moved to the external lay- 
ers of the HMNS via the magnetic-field tension. When this 
happens, the centrifugal support becomes insufficient to bal- 
ance the gravitational forces and the HMNS is induced to col- 
lapse to a rotating BH with dimensionless spin J/M 2 ~ 0.80 
(cf. Table [fl}, Such oscillations are directly related to the os- 
cillations visible in the evolution of the maximum of the rest- 
mass density in the top-left panel of Fig. [2] 



B. Low-mass binaries 

In analogy with what was done for the high-mass binaries, 
we show in Fig. [9] the GW signal for the low-mass models 
and also in this case the bottom right panel shows both the 
M1.45-B12 (black solid line terminated at t - r ps 19ms) 
and Ml. 4 5 -B0 (red dashed line) models for comparison. 
Since we have not evolved these models until the collapse of 
the HMNS to BH, only the inspiral and the post-merger phase 
(the part of the signal for t— r > 8.5 ms) are present in the GW 
signal. The high-frequency oscillations in the post-merger 
phase are related to the formation of a bar-deformed HMNS 
(as already described in J2]), whose spinning frequency is 
not significantly affected by the presence of magnetic fields. 
Also in this case, all the waveforms are very similar to each 
other both during the inspiral and after the merger. As a re- 
sult, and in contrast with what was seen for the high-mass 
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FIG. 8: Gravitational waves for the high-mass binaries as a function of the retarded time t — r in ms. The last panel shows for comparison 
also the unmagnetized model (i.e., red dashed line which collapses at t — r w 17ms) together with the model Ml . 62-B12 (black solid line 
which collapses earlier). 



case, the differences in the phase evolution are very small, 
at least over the timescales considered here (cf. bottom right 
panel). Clearly, if the HMNS continues to exist for longer 
times (on the radiation-reaction timescale), then the small dif- 
ferences may grow sufficiently and lead to a detectable differ- 
ence. While the numerical simulation of the secular evolution 
of the HMNS represents a challenge that we will address in 
future work, its impact on the detectability of the magnetic 
field will be further discussed in the next Section. 



C. Detectability of the magnetic field 

In order to assess the possibility of distinguishing between 
the different waveforms and hence establish whether differ- 
ent magnetizations of the HMNS can be measured, we have 
computed the power spectral densities of the GWs discussed 
before and plotted them in Fig. [10] against the sensitivity 
curves of different ground-based GW detectors. In particu- 
lar, we show the scaled power spectral densities ft-+(/)/ 1/ ' 2 
for the high-mass case (left panel) and low-mass case (right 
panel) with an initial magnetic field (dot-dashed red line) of 
B « 10 8 G (first row), B « 10 10 G (second row), and 
B m 10 12 G (third row). In all the panels the sources are 
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FIG. 9: Gravitational waves for the low-mass binaries as a function of the retarded time t — r in ms. The last panel shows for comparison also 
the unmagnetized model {i.e., red dashed line which terminates at t — r — 20 ms) together with the model 1 . 4 5-B12 (black solid line which 
terminates earlier). 



considered to be placed at a distance of 100 Mpc. We report 
also the spectral densities for binaries without magnetic field 
(solid black line) to aid in the comparison. Also shown are 
the noise curves of the Virgo detector (dotted green line), of 
the advLIGO O and adv Virgo J47) detectors (short-dashed 
blue and dot-dashed magenta lines, respectively), and of the 
planned Einstein Telescope [15 48 1 (dashed red line). With 
a dotted vertical line we indicate the value of twice the initial 
orbital frequency /q, so that the signal for / < 2/o should be 
ignored. 

In all the panels the part of the signal up to about 700 Hz 
is associated with the inspiral part of the waveform and in the 
case of the high-mass binaries (panels in the left column) it is 



also the strongest peak. The low-mass binaries (panels in the 
right column) also show an additional peak with an amplitude 
comparable to that at / 700 Hz and it is related to (twice) 
the spinning period of the bar-deformed HMNS. That peak 
appears for all the models at a frequency of approximately 
2 kHz and its amplitude is sufficiently high to enter into the 
band of advLIGO. 

While we expect the position in frequency of the peak to 
be accurate, its amplitude clearly depends on the subsequent 
evolution of the HMNS, which we have followed here only 
for about 12 ms. Clearly, should the HMNS survive on much 
longer timescales as shown in 1 8 1 (see the right panel of Fig. 
Al in the Appendix of [8|), then the amplitude of this peak 
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FIG. 10: Scaled power spectral densities h+ (f)f 1 ^ 2 for the high-mass case (left panel) and low-mass case (right panel) without magnetic field 
(solid black line) or with an initial magnetic field (dot-dashed red line) of « 10 8 G (first row), « 10 10 G (second row) or « 10 12 G (third 
row). In all the panels the sources are considered when placed at a distance of 100 Mpc. Shown also are the noise curves of the Virgo detector 
(dotted green line), of the advanced LIGO and advanced Virgo detectors (short-dashed blue and dot-dashed magenta lines, respectively), and 
of the planned Einstein Telescope (dashed red line). The dotted vertical lines indicate the value of twice the initial orbital frequency fo- 



16 



could be considerably larger and could scale with the square 
root of the period in which the HMNS continues to rotate be- 
fore collapsing. Even when leaving aside the role that the 
energy extraction via neutrinos may play on the evolution of 
the post-merger object, the hydrodynamical survival of the 
bar deformation in a rapidly rotating star is still a matter of 
debate. The general-relativistic simulations of isolated NSs 
first carried out in Ref. ||49l and then analyzed in great detail 
in Ref. ll50H52l . all indicate that the bar deformation persists 
only over a timescale which is comparable with (or slightly 
larger than) the dynamical one 6 . This is due to the coupling 
between the m = 2 bar deformation with other unstable 
modes (most notably the m = 1 one), which grow to com- 
parable amplitudes and suppress the instability, redistributing 
angular momentum (see also fl55ll for a perturbative analysis 
in terms of a Faraday resonance). On the other hand, sim- 
ulations of stellar-core collapse (see [56 1 for a recent review 
and a complete set of references) and the very long simula- 
tions carried out in J8J suggest that bar-deformed stellar cores 
or HMNSs can be produced and survive on timescales much 
longer than the dynamical one. This different behavior in the 
persistence of the bar deformation may well be due to the very 
different distribution of angular momentum and density strat- 
ification between the two configurations. Work is ongoing to 
confirm whether this is actually the case. 

Also quite evident from all the panels is that the spectra are 
very similar but not identical and that these differences be- 
come more appreciable for larger initial magnetic fields. In- 
deed the largest differences appear for B ~ 10 10 and, as for 
the accelerated collapse discussed in Fig. [5] magnetic fields 
of this strength are those that most influence the postmerger 
dynamics. Once again, it is worth emphasizing that the spec- 
tra presented here refer to a possibly too short portion of the 
evolution of the HMNS and if the HMNS does survive on 
much longer timescales, then the small differences shown here 
would become considerably more pronounced and well within 
the sensitivities of advanced detectors. 

In order to asses in a more quantitative way the possibility 
to detect these small differences in the GWs, we have com- 
puted the overlap between two waveforms h B1 , h B2 from bi- 
naries with initial magnetic fields Bl, B2 as 



0[h B1 ,h B 



(KAK 



where (h B1 \h B2 ) is the scalar product, defined as 



(16) 



df 



s h (f) 



(17) 



and h(f) is the Fourier transform of the GW h(t) and Sh(f) 
is the noise power spectral density of the detector (we have 
considered advLIGO here). Taking two waveforms, the closer 



6 Similar results have been found also in Newtonian simulations 1531 and 
also for magnetized stars 1541 



their overlap is to 1, the harder will be for a detector to distin- 
guish them. 

The overlaps computed for all the magnetized binaries con- 
sidered here when compared with the corresponding non- 



magnetized models are collected in Table III Note that we 
present both the total overlap O, i.e., the overlap computed 
over the full time-series, and the overlaps computed over the 
inspiral only or the post-merger only, i.e., mS p and O 



respectively. Given the values in Table III and since present 



postm? 



and advanced detectors could potentially distinguish two sig- 
nals if O < 0.995, it is clear that a detector such as advLIGO 
or adv Virgo would not be able to distinguish between a mag- 
netized binary and an unmagnetized one (cf. second column in 
the Table). Similar considerations apply also when the over- 
lap is computed only over the inspiral phase (cf. third column 
in the Table). However, if the overlap is computed only over 
the post-merger phase (cf. fourth column in the Table) then 
it is evident that the differences among the various binaries 
are much larger and the corresponding overlaps considerably 
smaller. Hence, we conclude that a long-lived HMNS and a 
detector with sufficient sensitivity at high frequencies (such as 
the Einstein Telescope) could be able to measure the level of 
magnetization in the progenitor NSs. 

To complete the information about the GW emission from 
magnetized BNSs, we have also computed the signal-to- 
noise-ratio (SNR) defined as 



1M/)1 ; 
s h (f) 



df 



(18) 



for different detectors and we have listed their values in Ta- 



ble III for a source at 100 Mpc. Overall, it is easy to realize 



that while the current Virgo and LIGO detectors (respectively 
the fifth and sixth columns) would not be able to detect these 
signals, SNRs larger than 1 are obtained when considering ad- 
vLIGO and advVirgo, and even larger than 40 in the case of 
the Einstein Telescope (last column in the Table). It is worth 
stressing that these SNRs should be seen as lower limits. First, 
the binaries are expected to enter the sensitivity band at lower 
frequencies than the ones considered here, hence adding con- 
siderable power to the SNR. Second, as discussed extensively 
above, the possibility of a long-lived HMNS could signifi- 
cantly add to the power at high frequencies, hence increasing 
the SNR. 

In summary, the results presented here indicate that BNSs 
do represent strong sources of GWs and that these can be de- 
tected at distances up to 100 Mpc by the planned advanced 
interferometers. Determining the level of magnetization of 
the progenitor stars will be very difficult if the detected sig- 
nal is confined essentially to the inspiral, while it could be 
possible if the HMNS survives for sufficiently long times as a 
deformed and spinning bar. In this latter case, detectors which 
have high sensitivities at high frequencies, such as advLIGO 
and more importantly the Einstein Telescope, will be in a good 
position to measure the strength of the magnetic fields and 
hence extract important physical and astrophysical informa- 
tion on the progenitor NSs. 
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TABLE III: GW-related quantities. Column 2 shows the total overlap computed for advLIGO between the magnetized models and the cor- 
responding unmagnetized binary, while columns 3 and 4 represent the overlap computed over the inspiral and over the post-merger phase. 
Finally, columns 5 — 9 show the SNR computed for different detectors for all the eight models considered here. The SNR has been obtained 
assuming a source at 100 Mpc. 



Binary 




o 


C^insp 


C^postm 


SNR (Virgo) 


SNR (LIGO) 


SNR (adv Virgo) 


SNR (advLIGO) 


SNR (ET) 


Ml . 45- 


B0 


1.000 


1.000 


1.000 


0.33 


0.23 


1.94 


2.11 


38.90 


M1.45- 


B8 


0.997 


0.999 


0.926 


0.33 


0.23 


1.94 


2.10 


38.72 


M1.45- 


BIO 


0.996 


0.999 


0.934 


0.33 


0.23 


1.94 


2.11 


38.82 


M1.45- 


B12 


0.996 


0.999 


0.899 


0.33 


0.23 


1.94 


2.11 


39.01 


Ml . 62- 


BO 


1.000 


1.000 


1.000 


0.36 


0.25 


2.00 


2.24 


42.57 


Ml . 62- 


B8 


0.998 


1.000 


0.938 


0.36 


0.25 


2.00 


2.24 


42.59 


Ml . 62- 


BIO 


0.993 


1.000 


0.724 


0.36 


0.25 


2.00 


2.23 


42.48 


Ml . 62- 


B12 


0.997 


1.000 


0.893 


0.36 


0.25 


2.00 


2.24 


42.49 



V. CONCLUSIONS 

There is little doubt that BNSs represent prime sources for 
present and advanced GW detectors. Equally clear is that 
NSs are observed to have large magnetic fields, with values 
which can be as high as 10 16 G for isolated and young magne- 
tars. It is therefore of great importance to assess what role the 
magnetic fields play during the inspiral and merger on BNSs. 
Extending the research presented in 0, we have presented 
the first numerical simulations of magnetized BNSs with as- 
trophysically realistic magnetic fields. More specifically, we 
have carried out a systematic investigation of the dynamics of 
both matter and magnetic fields of equal-mass BNSs. While 
previous works ||3]|5]|6) considered only astrophysically unre- 
alistic magnetic fields (B R3 10 16 — 10 17 G) or focused mainly 
on the inspiral part [6j, here we have considered magnetic- 
field values ranging from 10 s to 10 12 G, and evolved BNSs 
through all the stages of the inspiral, merger, HMNS evolu- 
tion, and collapse to BH. 

Overall, we have shown that realistic magnetic fields do not 
affect sensibly the dynamics of the inspiral, but they can in- 
fluence that of the post-merger, where they can accelerate the 
collapse of the HMNS. The different time intervals from the 
merger to the collapse of the HMNS also imply that the tori 
produced around the BH have slightly different masses, re- 
flecting the different distributions of matter and angular mo- 
mentum at the time of collapse. As a result of the tight corre- 
lation between the degree of magnetization of the NS matter 
and the delay time of the collapse, the measurement of the lat- 
ter via a GW detection will allow us to infer the former. To 
the best of our knowledge, this is the first time that effects of 
this type have been discussed in the evolution of inspiralling 
and magnetized NSs. 

Magnetic fields can be amplified at the merger of the bi- 
nary, when a Kelvin-Helmholtz instability develops between 
the outer layers of the two stars. Although the resolution used 
here is the highest employed so far in simulating magnetized 
BNSs and it is sufficient to reveal the development of the in- 
stability and the exponential growth of the toroidal magnetic 
field, the amplification we have measured is only of about one 
order of magnitude and is much smaller than that reported in 



Ref. l40l . where the newly produced fields reach values in 
equipartition with the kinetic energy. Although it is possi- 
ble that the different results are due to the different numerical 
methods employed in Ref. [40 1, we believe the reason behind 
our modest amplifications to be that the shear layer between 
the two stars survives only for about 1 ms, before being de- 
stroyed by the collision between the two stellar cores. Such a 
short timescale and the relatively small velocities at the shear 
layer are probably insufficient to yield the type of amplifica- 
tion that has been obtained in more "controlled" simulations 
of the Kelvin-Helmholtz instability BT1I421 . 

The toroidal magnetic field continues to be amplified also 
after the Kelvin-Helmholtz instability has been suppressed 
and it can reach values that are comparable with the initial 
poloidal one either during the evolution of the HMNS (in the 
case of low -mass binaries) or during the evolution of the torus 
produced after the HMNS collapse (in the case of high-mass 
binaries). This result is particularly important since it suggests 
that the magnetic -field topology in the tori formed from BNS 
mergers is not purely poloidal, contrarily to what has been 
assumed so far by some other research groups that perform 
simulations of magnetized accretion disks. 

When considered in terms of their GW emission, the mag- 
netized binaries studied here show that it is unlikely that the 
degree of magnetization will be measurable by present and 
advanced detectors if the inspiral is the only part of the sig- 
nal available. However, if the HMNS survives for sufficiently 
long times as a deformed and spinning bar, then the modi- 
fications introduced by the presence of magnetic fields could 
lead to waveforms which differ appreciably from those of non- 
magnetized binaries. In this case, detectors which have high 
sensitivities at frequencies larger than about 2 kHz, such as 
advLIGO and, more importantly, the Einstein Telescope, will 
be able to measure these effects for binaries up to distances of 
about 100 Mpc. 
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